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ABSTRACT 

Due  to  the  lack  of  longitudinal  stiffening  along 
submarine  hulls,  they  are  inherently  vulnerable  to 
underwater  explosions  directly  off  the  bow  or  stern.  ' 
Accordion-like  deformations  of  the  hull  are  set  into  motion 
which  could  cause  dynamic  amplification  in  the  transient 
response  of  attached  substructures.   In  underwater  shock 
acceptance  tests  of  internal  equipment,  this  interaction  is 
created  by  exploding  a  charge  in  a  fore  and  aft  configura- 
tion with  the  submerged  shock  test  vehicle  (SSTV) .   With  the 
increasing  availability  of  large  computers  and  the  rapid 
development  of  numerical  methods,  several  computer  codes 
have  been  written  to  predict  equipment  response  to 
underwater  shocks.   Using  the  ELSHOK  (ELASTIC  SHOCK)  code, 
this  investigation  studies  the  effect  of  hull/substructure 
interaction  on  stiffened  shell  response  at  resonance 
following  an  end-on  load.   The  transient  response  of  the  • 
coupled  shell/substructure  system  from  tapered  and 
conventional  charges  of  equivalent  impulse  is  examined  in 
this  study. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Submarines  are  inherently  vulnerable  to  explosions 
detonated  directly  off  the  bow  or  stern.   While  they  have 
adequate  protection  against  side  loads  with  transverse 
bulkheads  and  frames,  there  is  little  longitudinal 
protection  along  the  hull.   As  a  result  of  this,  they  are 
susceptible  to  accordion-like  deformations  of  the  hull  upon 
impact  from  end-on  loads.   It  is  possible  that  the  frequency 
of  the  hull  motions  could  excite  internal  equipment  into 
resonance  thereby  causing  unacceptable  damage.   This 
investigation  examines  the  coupled  elastic  response  of  a 
shell/substructure  system  when  subjected  to  an  end-on 
underwater  explosion  (UNDEX)  as  the  mass  and  stiffness  of 
the  internal  substructure  is  varied. 

The  current  specifications  for  conducting  underwater 
shock  tests  on  submarine- installed  equipment  are  contained 
in  MIL-S-901D  [Ref.  1].   This  document  specifies  the 
explosive  charge  weight  and  geometry  of  the  test,  as  well  as 
the  mounting  and  orientation  of  the  equipment  being  tested 
within  the  submerged  shock  test  vehicle  (SSTV) .   The  SSTV  is 
nothing  more  than  a  ring-stiffened  cylindrical  shell  with 
circular  endplates  designed  to  simulate  the  hull  motions 
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created  by  an  underwater  shock.   Recognizing  the  potential 
for  unacceptable  damage  to  hull  mounted  equipment  from  an 
end-on  load,  the  specification  requires  the  first  of  four 
UNDEX's  to  be  performed  in  a  fore  and  aft  configuration. 
The  suitability  of  an  equipment  design  or  installation  is 
evaluated  according  to  its  ability  to  function  as  intended 
during  and  after  each  shock  impulse.   Equipment  tests  of 
this  nature  are  quite  expensive  and  require  a  great  deal  of 
preparation.   It  is  desirable  for  the  designer  to  have  some 
idea  of  what  the  transient  response  of  the  hull  and 
equipment  will  be  prior  to  the  actual  UNDEX. 

The  development  of  analytical  methods  and  computer  codes 
for  the  analysis  of  the  reaction  of  submerged  structures  to 
underwater  explosions  has  enabled  the  engineer  to  simulate 
the  shock  response  of  submarines  and  equipment  with 
increasing  accuracy.   Fini te-element/'f ini te  difference 
methods  allow  for  the  analysis  of  structure  and -fluid 
responses,  and  the  work  of  Geers  [Ref.  2]  provides  a  means 
to  incorporate  fluid-structure  interaction  effects.   Several 
computer  codes  currently  exist  which  utilize  these 
principles  to  analyze  installed  equipment  response  to  shock 
waves.   One  of  them  is  the  ELSHOK  (ELASTIC  SHOCK)  code 
developed  by  Weidlinger  Associates,  under  the  sponsorship  of 
the  Defense  Nuclear  Agency  (DNA)  and  the  Office  of  Naval 
Research  (ONR) .   This  code  was  developed  to  investigate 
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modern  submarine  underwater  explosive  shock  response  in 
conjunction  with  a  testing  program  using  small  to  large 
scale  models  and  shaped  or  tapered  explosive  charges.   The 
accuracy  of  ELSHOK  has  been  validated  in  several  highly 
controlled  tests,  notably  the  1983  low  level  explosive  test 
of  an  SSN  668  class  submarine.   The  code  was  made  available 
to  the  Naval  Postgraduate  School  by  the  DNA  with  support 
from  Weidlinger  Associates.   It  was  first  used  at  the  school 
by  LT  Mark  Welch,  USN,  in  his  study  comparing  the  shock 
response  predictions  of  ELSHOK  with  those  obtained  using  the 
Dynamic  Design  Analysis  Method  (DDAM)  [Ref.  3]. 

B.   PURPOSE  FOR  THIS  INVESTIGATION 

An  explosive  detonated  underwater  exerts  a  great  amount 
of  pressure  on  the  surrounding  fluid.   The  water  is  actually 
compressed  by  the  force  of  the  explosion  causing  a  shock 
wave  to  form  at  the  point  of  detonation  which  propagates  in 
a  roughly  spherical  shape  at  the  speed  of  sound  in  water. 
When  the  pressure  pulse  impinges  on  a  submarine,  dynamic 
responses  result  in  the  hull  from  the  fluid-structure 
interaction.   The  submarine  experiences  translat ional  motion 
away  from  the  point  of  impact.   The  velocity  of  this  motion 
is  of  particular  interest  to  a  designer  of  submarines  or 
weapons . 

In  an  explosion  with  a  side-on  aspect  with  the 
submarine,  the  shock  wave  hits  the  hull  in  the  middle  before 
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it  hits  the  ends  due  to  its  spherical  shape.   This  causes 
the  middle  to  bow  away  from  the  point  of  impact  while  the 
ends  remain  initially  fixed.   The  resulting  motion  is  called 
"whipping"  as  the  hull  bends  back  and  forth  about  the 
middle.   With  an  explosion  off  the  bow  or  stern,  the  shock 
wave  hits  the  submarine  as  a  plane  wave  due  to  the  small 
cross-sectional  area  at  the  point  of  impact.   Since  the 
force  of  the  shock  wave  is  axisymmetric  about  the  centerline 
of  the  submarine,  no  whipping  motion  is  induced  in  the  hull. 
The  inertial  forces  are  transmitted  down  the  hull  until  they 
reach  the  stern.   The  stern  reacts  violently  to  the  combined 
inertial  forces  resulting  from-  the  explosion  as  well  as  the 
reflected  pressure  pulses  from  the  water  displaced  by  the 
translat ional  motion  of  the  submarine.   After  the  initial 
transient  response,  the  hull  settles  into  an  "accordion 
mode"  about  the  middle  where  the  motion  of  the  bow  opposes 
that  of  the  stern. 

The  purpose  of  this  investigation  is  to  examine  the 
shell/substructure  interaction  when  subjected  to  UNDEX  end- 
on  loading  using  the  ELSHOK  code.   The  hull  used  in  this 
study  is  a  ring-stiffened  cylindrical  shell  with  endplates 
similar  in  shape  to  those  used  as  SSTV's.   The  internal 
equipment  or  substructure  is  a  diaphragm  attached  to  the 
shell  whose  thickness  is  increased  in  order  to  study  how  the 
mass  ratio  between  substructure  and  shell  affects  the 
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coupled  transient  response.   The  dynamic  amplification  of 
the  response  of  the  substructure  at  resonance  is  examined 
for  both  taper  and  conventional  charges  of  equivalent 
impulse. 


16 


II.   EXPLOSIVE  SHOCK  RESPONSE  USING  ELSHOK 

A.   GENERAL  PRINCIPLES  OF  OPERATION 

The  ELSHOK  computer  code  consists  of  a  family  of 
programs  developed  to  calculate  the  transient  response  of  a 
submerged,  ring-stiffened  shell  of  revolution  of  finite 
length  to  an  underwater  shock  wave  emanating  from  an 
explosive  source  located  at  an  arbitrary  point  away  from  the 
structure  [Ref.  4].   The  shell  is  assumed  to  be  linearly 
elastic,  with  or  without  internal  substructures,  and  the 
surrounding  fluid  is  treated  as  an  infinite  acoustic  medium. 
Component  modal  analysis  is  employed  in  all  phases  of  the 
calculations.   The  complete  structural  system  is  considered 
to  consist  of  the  ring-stiffened  shell  and  any  attached 
substructures.   The  vibration  modes  of  each  component  are 
calculated  sepa-rately,  and  the  equations  of  motion  for  the 
entire  system  are  obtained  by  enforcing  compatibility  of 
deformation  at  the  points  of  attachment.   The  free-free 
modes  of  the  empty  ring-stiffened  shell  and  the  fixed-base 
modes  of  each  individual  substructure  are  coupled  through 
the  use  of  dynamic  boundary  conditions  [Refs.  5,  6].   This 
eliminates  the  need  to  calculate  the  combined  modes  and 
natural  frequencies  of  the  entire  system  as  well  as  the 
requirement  for  a  combined  system  stiffness  matrix. 
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The  structure-fluid  interaction  is  approximated  in 
ELSHOK  using  the  Doubly  Asymptotic  Approximation  (DAA) 
methods  of  Geers  [Ref.  2],    expressed  in  terms  of  functions 
which  are  orthogonal  over  the  wet  surface  of  the  submerged 
shell.   By  matching  exact  pressure-velocity  relations  at 
zero  and  infinite  frequencies,  the  elements  of  the  matrices 
in  the  resulting  DAA  are  obtained.   In  transient  problems, 
the  DAA  yields  exact  solutions  at  early  and  late  times 
providing  a  smooth  transition  between  these  two  limits.   By 
accounting  for  the  effects  of  the  fluid  with  quantities 
defined  solely  on  the  wet  surface  of  the  shell,  the  DAA 
essentially  uncouples  the  fluid  field  from  the  structural 
field. 

B.   ORGANIZATION  AND  OPERATION  OF  ELSHOK 

As  stated  previously,  the  ELSHOK  code  contains  a  group 
of  computer  programs  which  are  utilized  in  order  to  obtain 
velocity-time  histories  for  various  locations  in  the  shell- 
substructure  system.   The  programs,  in  order  of  execution, 
are : 

1)  B0S0R4 — structural  analyzer  for  shell  [Ref.  7] 

2)  ACESNID--virtual  mass  processor 

3)  PIFLASH — shell-fluid  processor 

4)  SAPIV — structural  analyzer  for  substructure  [Ref.  8] 

5)  PICRUST--substructure  processor 

6)  USLOB--time  integration  processor 
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7)   PUSLOB — plotting  processor 
Figure  1  illustrates  the  general  relationship  between  the 
programs  which  make  up  the  ELSHOK  computer  code.   As  can  be 
seen,  the  operation  is  divided  into  four  phases. 
1.   Phase  I--Shell  and  Fluid  Analysis 

The  first  step  of  the  ELSHOK  calculation  is  to 
determine  the  in-vacuo  free-free  modes  and  natural  frequen- 
cies of  the  shell.   This  is  accomplished  by  modeling  the 
shell  using  the  B0S0R4  finite  difference  code.   This  code 
applies  to  segmented,  ring-stiffened,  branched  shells  of 
revolution  having  various  meridional  geometries,  wall 
constructions,  and  ring  reinforcements  making  it  well  suited 
for  the  modeling  of  submarines  or  SSTV's.   For  compatibility 
with  ELSHOK,  a  separate  B0S0R4  calculation  must  be  performed 
for  each  circumferential  harmonic  (N)  included  in  the 
analysis.   Examples  of  different  harmonics  in  the  circum- 
ferential distributions  include  M  =  0  (breathing/torsional) 
and  N  =  1  (pure  translation/whipping)  as  shown  in  Figure  2. 

The  second  part  of  Phase  I  is  the  computation  of  the 
virtual  mass  array  using  ACESNID.   This  provides  the 
late-time  contribution  of  the  DAA.   The  virtual  mass  array 
is  determined  from  the  solution,  based  on  simple  sources  of 
a  low-frequency  steady-state  problem  in  which  normal 
displacements  corresponding  to  surface  expansion  functions 
are  appli-ed  to  the  surface  of  revolution  in  the  infinite 
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Figure  2.   Nodal  Patterns  for  Cylindrical  Shells 
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fluid  having  the  same  shape  and  size  as  the  net  surface  of 
the  shell.   ACESNID  requires  only  one  execution  for  all 
values  of  N  considered. 

The  execution  of  PIFLASH  completes  Phase  I  by 
combining  the  outputs  from  B0S0R4  and  ACESNID  to  create  a 
"shell-fluid  file"  which  contains  all  the  information 
required  to  describe  the  shell  during  subsequent  operations. 
One  of  the  useful  features  of  ELSHOK  is  the  fact  that  any 
number  of  substructures  can  be  analyzed  for  a  given  shell 
without  having  to  repeat  Phase  I  calculations. 

2.  Phase  II--Substructure  Analysis 

The  analysis  of  the  substructure  is  performed  using 
the  SAPIV  finite  element  code.   The  SAPIV  code  contains  a 
variety  of  finite  elements  (plates,  beams,  pipes,  etc.) 
which  can  be  used  in  the  modeling  process.   It  is  a  general 
purpose  code  based  upon  a  lumped  mass  formulation  applying 
to  linearly  elastic  structures.   The  PICRUST  code  takes  data 
from  SAPIV  and  reorganizes  it  to  facilitate  the  solution  of 
the  equations  governing  the  transient  response  problem  under 
study.   It  is  during  the  execution  of  PICRUST  that  the 
connectivity  between  the  substructure  and  the  shell  is 
accounted  for  and  influence  coefficients  for  evaluating  the 
forces  developed  at  the  points  of  attachment  are  calculated. 

3 .  Phase  1 1 1--Submerged  Shock  Response 

In  order  to  prepare  a  combined  shell-substructure 
input  for  the  USLOB  time  integration  processor,  the 

22 


shell-fluid  file  from  Phase  I  is  merged  with  the 
substructure  files  from  Phase  II.   USLOB  is  the  portion  of 
the  ELSHOK  code  which  contains  the  underwater  shock  from  an 
arbitrary  point.   It  allows  pressure-time  history  inputs  for 
taper  charge  modeling  or  an  exponentially  decaying  pressure 
impulse  such  as  those  encountered  with  conventional  charges. 
It  is  of  the  form  [Ref .  9] : 

P(t)  =  K^{^^^^/R)     ^    expi-t/Q^)  (1) 

where : 

P(t)  =  incident  pressure  on  the  shell  (psi) 

K,  =  multiplicative  constant  for  incident  pressure 

K2  =  spatial  decay  constant  for  incident  pressure 

t     =  time  after  arrival  of  shock  wave  at  point  of 
interest  (msec) 

"w     =  weight  of  spherical  charge  (lb.) 


'o 


1/3   1/3     4 
=  K^w  ^     {w   ^    /R)     =  time  constant  of  exponential 

decay  (msec) 


K^    =  multiplicative  constant  for  time  constant 

K.    =  spatial  decay  constant  for  time  constant 

R     =  distance  from  the  explosive  to  the  point  of 
interest  (ft) 

K, ,  K2f  K^,  and  K^  are  constants  which  depend  on 

explosive  type.   USLOB  employs  a  modified  version  of  the 

Runge-Kutta  integration  method  to  produce  velocity-time 

histories  in  tabular  form  at  user  specified  points  on  the 

shell  and  the  substructure. 
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4 .   Phase  IV--Plots  of  Velocity-Time  Histories 
After  the  velocity-time  histories  have  been 
obtained,  PUSLOB  is  utilized  to  display  the  velocity 
responses  in  a  plotted  format  on  a  TEKTRONIX  graphics 
terminal.   These  plots  are  adequate  for  viewing  trends  and 
general  responses,  however  for  formal  display  purposes,  the 
IBM  Versatec  plotter  yields  a  superior  product.   Appendix  A 
contains  a  code  which  was  written  to  convert  the  velocity 
punch-card  files  created  by  PUSLOB  into  a  data  format  which 
is  acceptable  for  Versatic  plotting  using  the  EASYPLOT 
program. 
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III.   MODEL  USED  IN  THE  ANALYSIS 

A.   SHELL  MODEL 

The  shell  used  in  this  analysis  is  a  small-scale  model 
of  a  typical  SSTV  used  for  conducting  underwater  shock  tests 
on  high-impact  shipboard  machinery.   As  can  be  seen  in 
Figure  3,  it  is  a  high  strength  steel  free-free  ring 
stiffened  cylinder  with  aluminum  endplates.   Appendix  B 
includes  a  sample  B0S0R4  input  for  one  of  the 
circumferential  harmonic  distributions.   Data  can  be  input 
to  B0S0R4  using  a  formatted  input  code  or  by  responding 
interactively  to  a  set  of  prompts.   Each  segment  is  modeled 
separately  which  allows  for  different  materials  and  physical 
properties  to  be  included  in  the  structure.   The  small 
stiffeners  along  the  shell  are  represented  by  an  orthotropic 
approximation  which  in  effect  increases  the  density  of  the 
shell.   The  six  discrete  rings  are  modeled  individually, 
allowing  for  different  size  or  material.   The  torsional 
rigidity  (GJ)  for  each  ring  is  found  using  the  following 
relation : 

GJ  =  EJ/2  (1  +  V).  (2) 

where  : 

E  =  modulus  of  elasticity  (psi) 
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4 
J  =  A  /(40  I  )  =  torsional  rigidity  constant  [Ref.  10] 

ir 

V  =  Poisson's  ratio 

Since  the  shell  is  a  surface  of  revolution,  the 
properties  of  each  ring  cross-section  need  only  be  specified 
at  one  point.   The  end  plates  (segments  one  and  three)  are 
divided  into  ten  equally  spaced  nodes,  and  the  cylindrical 
portion  of  the  shell  (segment  two)  is  divided  into  forty- 
five  nodes  along  the  longitudinal  axis.   As  with  any  finite 
difference  code,  greater  accuracy  will  be  obtained  in  the 
solution  with  higher  numbers  of  nodes,  however,  computa- 
tional time  will  increase.   Two  additional  nodes  are 
automatically  inserted  by  B0S0R4  into  each  segment  in  order 
to  reduce  the  truncation  errors  associated  with  segment 
interfaces  and  to  prevent  spurious  vibration  or  buckling 
modes.   Through  finite  difference  techniques,  B0S0R4 
determines  the  in-vacuo  free-free  modes  and  natural 
frequencies  of  the  shell.   Since  this  study  deals  strictly 
with  end-on  loading,  the  problem  is  axisymmetric  about  the 
longitudinal  axis.   Although  only  the  N  =  0  (breathing/  • 
torsional)  modes  are  activated,  a  few  N  =  1  ( translat ional ) 
models  must  be  calculated  and  retained  due  to  the  fact  that 
the  ELSHOK  computer  code  was  developed  for  side-on  loading 
and,  therefore,  requires  a  translat ional  input.   Since  there 
is  no  translation,  the  solution  is  not  affected  by  these 
superfluous  modes.   B0S0R4  generates  a  shell  file  which 
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gives  the  mode  shapes  for  each  natural  frequency.   All 
breathing  modes  have  been  retained  for  further  study, 
however,  the  torsional  modes  have  been  dropped  from  the 
analysis  since  they  are  not  activated  by  an  end-on  load. 
Table  I  contains  a  listing  of  modes  and  natural  frequencies 
retained  for  input  into  the  PIFLASH  program. 

Figure  4  gives  a  representation  of  a  typical  shell- 
substructure  configuration.   The  two  coordinate  systems  used 
in  an  ELSHOK  analysis  are  illustrated.   X,  Y,  and  Z  refer  to 
the  global  coordinate  system  of  the  shell,  while  x,  y,  and  z 
refer  to  the  local  coordinate  system  of  the  substructure. 
In  both  systems,  the  x-axes  run  longitudinally  down  the 
hull.   The  substructure  z-axis  is  related  to  the  global 
Z-axis  by  an  angle  a   where  a  refers  to  the  substructure. 
Points  on  the  wet  surface  of  the  shell  are  located  by 
specifying  the  meridional  arc  length  s,  measured  along  the 
reference  surface,  and  the  circumferential  coordinate  angle 
e.   The  lower  case  u,  v,  and  w  are  used  for  local  shell 
displacements.   To  eliminate  torsional  modes  from 
consideration,  all  N  =  0  modes  in  the  shell  file  with 
significant  v-d isplacements  are  discarded. 

Once  the  shell  files  have  been  obtained  from  B0S0R4,  the 
input  code  for  ACESNID  can  be  created.   This  program 
determines  the  virtual  mass  array  which  produces  the  late- 
time  contribution  of  the  DAA.   It  considers  a  cavity  in  an 
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TABLE  I 
MODES  AND  NATURAL  FREQUENCIES  OF  SHELL 

MODE  N  FREQUENCY  (HZ) 

2  0  1.259E-03 

4  0  2.239E+02 

6  0  3.850E+02 

7  0  4.761E+02 
9  0  5.935E+02 

12  0  7.826E+02 

13  0  9.369E+02 

14  _  0             _  1.048E+03 

19  0  1.584E+03 

20  0  1.585E+03 

21  0  1.596E+03 

22  0  1.598E+03" 

23  0  1.656E+03 

25  0  1.791E+03 

26  0  1.818E+03 

27  0  1.856E+03 

28  0  1.878E+03 

29  0  1.890E+03  • 

30  0  1.898E+03 

1  1  1.390E-04 

2  1  6.321E+01 
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infinite  acoustic  fluid  corresponding  to  the  wet  surface  of 
the  shell.   A  user  specified  number  of  surface  expansion 
functions  having  the  property  of  orthogonality  over  the  wet 
surface  of  the  shell  are  obtained  during  the  calculation  of 
the  virtual  or  entrained  mass.   The  cavity  is  divided  into 
bands  which  describe  the  behavior  of  the  fluid.   A 
sufficient  number  of  bands  must  be  provided  to  enable  the 
fluid  response  to  match  the  normal  motion  specified  by  the 
surface  expansion  functions.   In  a  sense,  they  can  be  looked 
upon  as  nodes  in  a  finite  element  model.   In  this  analysis, 

the  shell  is  assumed  to  be  totally  immersed  in  salt  water 

2    4 
with  a  mass  density  of  9.59684E-05  Ibf-s  /in  .   The 

endplates  are  divided  into  fifteen  bands  each,  and  the 

cylinder  has  fifty-one  bands.   Six  surface  expansion 

functions  have  been  generated  for  the  endplates  with  ten 

along  the  cylinder. 

The  shell  mode  files  from  B0S0R4  and  the  virtual  mass 

file  from  ACESNID  are  combined  in  the  PIFLASH  program  to 

produce  a  shell-fluid  file  for  further  processing.   The  • 

model  of  the  shell  is  now  complete  and  can  be  used  for  any 

variety  of  substructures  without  any  recalculation. 

B.   SUBSTRUCTUPE  MODEL 

The  idealized  internal  equipment  selected  for  study  in 
this  analysis  is  a  high  strength  steel  diaphragm  or  plate 
located  at  the  discrete  ring  at  frame  nine.   The  SAPIV 
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finite  element  code  is  utilized  to  model  the  diaphragm.   The 
diaphragm  has  been  divided  into  eight  "pie-shaped"  wedges 
which  are  further  subdivided  into  three  sections.   Twenty- 
four  plate  elements  are  used  in  this  model  with  twenty-five 
nodal  points.   Since  this  study  is  concerned  with  end-on 
loading,  the  node  at  the  center  of  the  plate  is  constrained 
to  move  only  in  the  x-direction  with  no  rotations.   The 
outer  nodes  (eighteen  through  twenty-five)  are  rigidly 
attached  to  the  shell.   This  is  accomplished  through  a 
modification  of  the  general  purpose  SAPIV  code  by  the 
developers  of  ELSHOK.   During  the  parametric  study,  the 
diaphragm  thickness  is  uniformly  varied  from  .25  inches  to 
9.25  inches.   Figure  5  is  a  TEKTRONIX  representation  of  the 
finite  element  model  of  the  substructure.   Appendix  B 
includes  one  of  the  SAPIV  input  codes  used  during  the 
parametric  study. 
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IV.   ANALYSIS 

The  purpose  of  this  analysis  is  to  study  the  shell/ 
substructure  interaction  after  the  system  is  subjected  to  an 
end-on  load  while  varying  the  mass  and  stiffness  of  the 
internal  substructure.   Two  different  types  of  loading  are 
explored:   a  taper  charge  and  a  conventional  charge.   Taper 
charges  are  characterized  by  sustained  incident  pressure  on 
a  body  over  a  period  of  time.   They  are  used  in  underwater 
shock  tests  to  simulate  the  type  of  pressure  profile  which  a 
body  would  experience  from  a  nuclear  detonation.   Extensive 
use  of  scaling  is  involved  in  efforts  to  reproduce  the 
physical  behavior  of  an  object  in  response  to  a  nuclear 
charge  without  generating  the  extremely  high  pressures 
associated  with  actual  detonations.   Conventional  charges 
create  a  large  initial  pressure  pulse  followed  by  a  rapidly 
decaying  exponential.   ELSHOK  uses  inputs  of  charge  weight 
and  standoff  geometry  to  calculate  the  transient  velocities 
of  the  shell  and  substructure  in  response  to  these  loads. 

A.   ANALYSIS  PROCEDURE 

Once  the  shell  model  has  been  produced,  a  series  of 
SAPIV  runs  are  conducted  for  the  diaphragms  of  increasing 
thickness.   Each  run  generates  natural  frequencies  and 
accompanying  mode  shapes  for  each  model.   The  PICRUST 
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program  reduces  the  data  obtained  through  SAPIV  into  a 
format  which  can  be  combined  with  the  shell-fluid  file  for 
processing  with  the  USLOB  code.   The  substructure  is 
connected  to  the  shell  during  this  phase  of  the  ELSHOK 
analysis.   The  diaphragm  is  attached  to  shell  segment  two  at 
frame  nine  (node  eleven) . 

The  weight,  type,  and  location  of  the  charge  is 
specified  in  the  USLOB  input  code.   To  set  up  an  end-on  load 
referenced  to  "the  global  coordinate  system,  the  charges  are 
placed  in  the  negative  X-direction  at  a  distance  of  70  feet 
(840  inches)  .   This  distance  is  typical  of  stand-off 
distances  used  during  actual  tests.   For  the  taper  charge,  a 
typical  incident  pressure-time  history  shown  in  Figure  6  is 
used  to  simulate  the  shock  wave  loading.   In  order  to  make  a 
meaningful  comparison  between  the  structural  responses  from 
a  taper  charge  and  a  conventional  charge,  it  is  necessary 
that  they  have  equivalent  shock  wave  impulses.   In  the  field 
of  underwater  shock  analysis,  "impulse"  refers  to  the  time- 
integral  of  the  pressure  profile  [Ref.  9].   The  impulse  of 
unit  area  of  the  shock  wave  front  up  to  a  time  t  after  its 
arrival  is  given  by: 


I(t)  =   /   P(t)dt 
-^0 


(3) 


The  impulse  of  the  taper  charge  is  found  by  calculating  the 
area  under  the  pressure-time  history  in  Figure  6.   As  was 
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Figure  6«   Incident  Pressure-Time  History — Taper  Charge 
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stated  earlier,  a  conventional  charge  follows  a  pressure- 
time  history  given  by  equation  (1) : 


where : 


P(t)  =  K^(w-'/-'/R)  ^  exp(-t/9^) 


K^w  ^  (w  /  /R)    =  time  constant  of  exponential 


(1) 


decay  (msec) 


HBX-1  is  a  popular  explosive  used  in  underwater  shock  tests, 
so  it  was  decided  to  use  it  in  this  study  for  the 
conventional  charge.   K,,  K^,    K^ ,  and  K.  are  constants  which 
depend  on  the  type  of  explosive  used.   The  values  used  in 
this  study  are  [Ref.  11]: 


K^  =  3.8354314E+05 


K2  =  1.144 


K^  =  3.03131E-05 


K^  =  -.247 


By  matching  the  impulses  between  the  taper  charge  and  the 
conventional  charge,  a  charge  weight  of  352  lbs  of  HBX-1  is 
found  to  yield  an  equivalent  impulse.   The  peak  pressure 
generated  by  the  charge  is  1620  psi  which  creates  the 
incident  pressure-time  history  shown  in  Figure  7.   The 
conventional  charge  expends  most  of  its  impulse  in  the  first 
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millisecond,  while  the  taper  charge  takes  over  four  milli- 
seconds to  expend  an  equivalent  amount  of  impulse.   Because 
of  the  sharp  rise  in  pressure  experienced  from  a  conven- 
tional charge,  a  higher  transient  response  is  developed  in 
the  shell  than  with  a  taper  charge  of  equivalent  impulse. 

After  selecting  suitable  charges,  the  time  step  incre- 
ment and  integration  limits  are  specified  in  the  USLOB  input 
code.   Through  a  trial  and  error  process,  enough  time  steps 
are  chosen  to  identify  the  significant  interactions  between 
the  shell  and  substructure.   It  was  found  that  1600  time 
steps  for  an  80  msec  time  period  gave  an  adequate  illustra- 
tion of  the  system  response.   The  USLOB  program  enables  the 
user. to  examine  the  velocity-time  history  response  of 
various  nodes  in  tabular  form.   This  gives  some  indication 
of  whether  the  run  was  successful,  however,  the  plots  from 
PUSLOB  are  required  before  response  characteristics  can  be 
identified. 

To  obtain  a  baseline  velocity-time  history,  a  complete 
analysis  is  performed  on  the  submerged  empty  stiffened 
shell.   From  the  velocity-time  history  plots,  the  dominant 
frequency  from  the  computed  shock  response  is  observed. 
Provided  the  mass  of  the  diaphragm  remained  small  with 
respect  to  the  shell,  it  was  felt  that  this  provided  a  good 
estimation  of  the  excitation  frequency  for  the  diaphragm. 
The  dominant  frequency  is  187.80  Hz  for  the  taper  charge  and 
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195.44  Hz  for  the  conventional  charge.   A  series  of  SAPIV 
calculations  are  conducted  on  diaphragms  of  varying 
thickness  in  order  to  identify  the  fundamental  frequency  of 
each  substructure.   These  values  are  given  in  Table  II.   An 
analysis  is  performed  on  the  shell/substructure  system  for 
each  diaphragm  thickness  to  observe  the  coupled  velocity- 
time  history  response  as  the  substructure  mass  grows  in 
relation  to  that  of  the  shell. 
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TABLE  II 
FUNDAMENTAL  FREQUENCIES  OF  DIAPHRAGM  SUBSTRUCTURE 

Thickness  (In)  Frequency  (Hz) 

.25  43.66 

.50  87.31 

.75  130.97 

1.00  174.59 

1.05  183.35 

1.06  185.10 

1.07  186.85 
1.10  192.10 
1.15  200.85 
1.25  218.20 
2.25  392.95 
3.25  567.55 
4.25  742.14 
5.25  916.73 

6.25  1091.33 

7.25  1265.92 

8.25  1440.51 

9.25  1615.42 
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V.   RESULTS 

The  velocity  time-history  response  begins  at  time  zero 
when  the  incident  pressure  wave  from  the  underwater 
explosion  strikes  the  shell.   With  end-on  loading,  the 
forward  endplate  receives  the  initial  contact  with  velocity 
indicated  along  the  longitudinal  axis.   The  empty  shell 
response  is  analyzed  first,  followed  by  a  discussion  of  the 
coupled  effects  when  the  diaphragm  achieves  resonance. 

A.   VELOCITY-TIME  HISTORY  RESPONSE — EMPTY  SHELL 

The  empty  shell  velocity-time  history  responses  for  the 
conventional  and  taper  charges  are  contained  in  Figures  8 
through  20.   It  can  be  seen  that  the  conventional  charge 
causes  a  much  greater  transient  response  in  the  shell.   In 
both  cases,  the  velocity  is  rapidly  damped  by  interaction 
between  the  shell  and  the  infinite  fluid  medium.   It  is 
interesting  to  observe  the  differences  in  response  between 
the  various  locations  on  the  shell.   The  forward  endplate  is 
perpendicular  to  the  shock  wave,  receiving  the  full  force 
from  the  incident  pressure  pulse.   Inertial  forces  are 
generated  by  the  endplate  causing  severe  local  deformation 
in  the  nearest  nodes  on  the  cylinder  due  to  the  rigid 
coupling  between  segments.   The  motion  propagates  down  the 
shell  along  the  longitudinal  axis  until  it  reaches  the  aft 
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endplate  over  a  millisecond  after  initial  shock  wave 
contact.   At  first  the  aft  endplate  moves  in  the  negative 
X-direction  as  the  "ripple"  proceeds  down  the  structure. 
Then  it  experiences  a  violent  transient  response  as  it  is 
affected  not  only  by  the  inertial  forces  from  the  shell,  but 
by  the  reflected  pressure  pulses  from  the  water  it 
displaces.   After  the  initial  shock,  the  shell  settles  into 
an  "accordion"  effect,  expanding  and  contracting  as  the 
forward  and  aft  endplates  move  in  opposite  directions. 

B.   VELOCITY-TIME  HISTORY  RESPONSE — COUPLED  SHELL/ 
SUBSTRUCTURE  SYSTEM 

From  the  results  of  the  analyses  performed  on  the  coupled 
shell/substructure  system.  Tables  III  and  IV  were  created  in 
order  to  identify  when  the  diaphragm  achieved  resonance.   i^/n 
is  the  frequency  ratio  between  the  fundamental  frequency  of 
the  simply-supported  diaphragm  (Q)  and  the  dominant  frequency 
observed  in  the  velocity-time  history  response  of  the  sub- 
merged empty  stiffened  shell  (fi) .   M1/M2  is  the  mass  ratio 
between  the  substructure  (M,)  and  the  shell  with  endplates 
(M2) .   Figure  21  is  a  plot  of  the  velocity  of  the  diaphragm 
center  versus  the  frequency  ratio  Q/Q,       It  reveals  that  in 
both  the  taper  and  conventional  charge  responses,  the 
diaphragm  achieves  its  maximum  velocity  or  resonance  when  the 
fundamental  frequency  of  the  substructure  is  nearly  equal  to 
the  excitation  frequency  from  the  shell.   Resonance  occurs 
when  the  mass  ratio  is  approximately  .08.   Since 
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TABLE  III 

MAXIMUM  VELOCITY  RESPONSE  OF  DIAPHRAGM  SUBSTRUCTURE 
WITH  VARYING  THICKNESS  (TAPER  CHARGE) 


Thickness  (In) 


n/U 


M1/M2 


Velocity  (ft/sec) 


.25 
.50 
.75 
1.00 
1.05 
1.06 
1.07 
1.10 
1.15 
1.25 
2.25 
3.25 
4.25 
5.25 
6.25 
7.25 
8.25 
9.25 


.232 

.018 

.465 

.036 

.762 

.055 

.930 

.073 

.976 

.076 

.986 

.077 

.995 

.078 

1.023 

.080 

1.069 

.084 

1.162 

.091 

2.092 

.164 

3.022 

.237 

3.952 

.310 

4.881 

.383 

5.811 

.456 

6.741 

.528 

7.670 

.601 

8.602 

,674 

22.19 

16.64 

20.79 

34.81 

35.65 

35.78 

35.65 

33.58 

30.36 

26.24 

11.61 

10.60 

9.04 

8.70 

8.44 

8.15 

7.87 

7.62 
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TABLE  IV 

MAXIMUM  VELOCITY  RESPONSE  OF  DIAPHRAGM  SUBSTRUCTURE 
WITH  VARYING  THICKNESS  (CONVENTIONAL  CHARGE) 


Thickness  (In)  n/n  M1/M2       Velocity  (ft/sec) 

.25  .223  .    .018  40.42 

.50  .447  .036  28.26 

.75  .670  .055  45.29 

1.00  .893  .073  68.73 

1.05  .938  .076  73.84 

1.06  .947  .077  '  74.71 

1.07  .956  .078  74.22 
1.10  . .  .983  .080  69.69 
1.15  1.028  .084  62.90 
1.25.  1.116  .091  62.75 
2.25  2.011  .164  21.77 
3.25  2.904  .237  17.67 
4.25  3.797  .310  15.51 
5.25  4.691  •  .383  13.69 
6.25  5.584  .456  12.01 
7.25  6.477  .528  11.37 
8.25  7.371  .601  11.00 
9.25  8.265  .674  10.69 
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Ji  =  n  =  Vk1/M1"   =  ^K2/M2"  (4) 

where:   Kl  =  overall  stiffness  of  the  substructure 
K2  =  overall  stiffness  of  the  shell 

it  can  be  said  that: 

Kl/Ml  =  K2/M2  =  K1/.08  M2 

or 

K1/K2  =  .08 

The  coupled  shell  and  substructure  can  then  be  modeled  as  a 
two  degrees  of  freedom  system  as  shown  in  Figure  22  where  Ml 
and  Kl  represent  the  substructure  and  M2  and  K2  represent 
the  shell.   The  large  mass  (M2)  has  an  uncoupled  natural 


frequency  (fj)  which  is  equal  to  v^K2/M2'.   This  is 
equivalent  to  the  dominant  frequency  (JT)  of  the  empty  shell. 
Since  this  is  a  two  degrees  of  freedom  system,  it  will 
possess  two  natural  frequencies  for  the  coupled  system.   It 
has  been  found  in  coupled  systems  that  for  decreasing  mass 
ratio  (M1/M2)  and  for  increasing  stiffness  ratio  (K1/K2) , 
the  first  and  second  modal  frequencies  of  a  two  degrees  of 
freedom  system  approach  the  decoupled  frequencies  f,  and  f^ 
of  two  single  degrees  of  freedom  systems,  respectively  [Ref. 
12].   The  first  natural  frequency  is  bounded  by  £2,    and  f-j^ 
is  bounded  by  the  second  natural  frequency.   From  Figure  23, 
it  is  seen  that  with  M1/M2  =  K1/K2  =  .08,  the  dominant 
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frequency  of  the  shell  is  approximately  equal  to  the  first 
modal  frequency  (f^i)  of  the  system.   The  result  is  that  the 
shell/substructure  system  achieves  resonance  when  excited  at 
the  dominant  frequency  of  the  empty  shell. 

With  the  diaphragm  at  resonance,  Figures  24  and  25 
illustrate  that  the  transient  response  of  the  diaphragm  is 
much  more  pronounced  than  that  of  the  shell  due  to  the  large 
differences  between  the  two  masses  as  well  as  the  soft 
spring  between  them.   Initially  the  diaphragm  is  drawn  in 
the  negative  X-direction  by  inertial  forces  as  the  shell  is 
displaced  by  the  shock  wave.   With  the  increased  motion  of 
the  shell  at  frame  nine,  the  diaphragm  receives  more  and 
more  energy.   At  approximately  t  =  20  msec,  the  diaphragm  is 
at  its  maximum  velocity  while  the  shell  is  relatively  calm. 
At  this  point,  the  diaphragm  begins  to  transfer  energy  back 
to  the  shell.   Since  kinetic  energy  is  equal  to  1/2  MV " ,  the 
diaphragm  is  a  small  mass  with  high  velocity  while  the  much 
larger  mass  of  the  shell  moves  at  a  lower  velocity.   The 
energy  of  the  system  is  transferred  between  the  shell  and 
substructure  as  the  overall  motion  is  damped  by  interaction 
between  the  shell  and  the  fluid  medium.   The  transfer  of 
energy  between  the  shell  and  the  diaphragm  creates  an  effect 
known  as  the  "beating"  phenomenon.   This  effect  is  seen  in 
the  velocity  time  history  responses  for  the  case  of 
diaphragm  resonance  in  Figures  26  through  39.   The  effect  is 
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quite  pronounced  at  the  endplates,  because  they  act  as  "hard 
spots"  receiving  the  energy  transmitted  through  the  shell. 
The  cylindrical  portion  of  the  shell  acts  as  a  conduit  for 
the  energy  transfer.   Again,  the  nodes  on  the  shell  nearest 
the  endplate  experience  some  local  deformation  due  to  the 
rigid  coupling  between  segments. 

The  empty  shell  velocity-time  history  response  is 
compared  with  the  resonant  diaphragm  case  in  Figures  40 
through  45.   Of  particular  interest  is  the  fact  that  the 
velocity  profiles  are  virtually  identical  for  the  first 
fifteen  milliseconds.   During  this  period  of  time,  the  shell 
is  receiving  energy  from  the  shock  wave  and  transferring  a 
portion  of  it  to  the  diaphragm.   Because  of  the  differences 
in  mass  between  the  shell  and  substructure,  not  much  energy 
is  lost  from  the  shell.   The  peak  velocities  do  not  vary  by 
more  than  five  feet  per  second  during  this  time  period.   The 
difference  arises  at  t  =  35  msec  when  the  shell  receives 
energy  from  the  resonating  plate.   The  decaying  velocity 
response  is  disturbed  during  the  exchange  of  energy  which-  is 
observed  in  the  beating  effect. 

As  the  thickness  of  the  diaphragm  is  increased  beyond 
the  resonant  condition,  the  interactions  between  the  shell 
and  substructure  are  reduced  dramatically.   The  shell/ 
substructure  velocity-time  history  remains  within  10%  of  the 
empty  shell  response  until  the  diaphragm  thickness  is  2.25 
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inches  thick  or  M1/M2  =  .164.   By  the  time  it  is  9.25  inches 
thick,  it  represents  over  half  the  mass  of  the  shell,  and  it 
acts  as  if  it  was  an  integral  part  of  the  shell  slaved  to 
the  motions  of  the  forward  endplate. 
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VI .   CONCLUSIONS 

Upon  inspection  of  the  velocity-time  history  responses 
at  various  locations  along  the  shell,  one  can  see  why  the 
potential  for  damage  is  so  great  from  an  underwater 
explosion  off  the  bow  or  stern  of  a  submarine.   While  the 
cylindrical  portion  for  the  most  part  serves  to  transmit 
inertial  forces  with  little  velocity  response,  the 
extremities  of  the  hull  experience  violent  transient 
responses  which  may  be  passed  to  attached  internal  equipment 
located  in  those  areas.   The  inertial  influence  of  the  mass 
of  the  endplates  on  neighboring  nodes  of  the  shell  is 
apparent  in  both  the  empty  shell  and  resonant  cases.   The 
response  from  a  conventional  charge  is  seen  to  be  much  more 
pronounced  than  that  of  a  taper  charge  of  equivalent  impulse 
due  to  the  rapid  expenditure  of  energy  against  the  forward 
endplate . 

The  results  obtained  from  the  resonant  diaphragm  case 
support  the  theory  that  with  decreasing  mass  ratio  between 
substructure  and  shell,  the  first  and  second  modal 
frequencies  of  the  two  degrees  of  freedom  system  approach 
the  decoupled  frequencies  of  the  substructure  and  shell  in 
two  single  degree  of  freedom  systems.   This  allows  the 
motion  of  hull  to  drive  the  diaphragm  into  resonance  when 
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the  dominant  frequency  of  the  shell  motion  is  roughly  equal 
to  the  fundamental  frequency  of  the  substructure.   It  is 
possible  that  this  could  occur  during  an  underwater  shock 
test  of  a  small  piece  of  equipment  in  a  large  SSTV.   This 
might  result  in  the  failure  of  the  equipment  in  a  scenario 
which  would  not  occur  in  a  real  submarine  with  its  myriad  of 
internal  masses  and  substructures. 

The  large  difference  in  mass  between  the  substructure 
and  the  shell  results  in  a  tremendous  dynamic  amplification 
in  the  response  of  the  diaphragm  at  resonance.   As  energy  is 
transferred  between  the  two  bodies,  the  large  mass  of  the 
shell  causes  its  velocity  to  be  quite  small  compared  to  -that 
of  the  diaphragm.   As  is  seen  in  Figures  40  through  45,  the 
initial  resonant  case  velocity  response  is  virtually 
unchanged  from  the  empty  shell  response.   If  the  mass  of  the 
internal  equipment  was  sufficiently  small  in  comparison  to 
an  SSTV  used  in  a  series  of  shock  tests,  one  ELSHOK  calcu- 
lation could  be  performed  on  the  empty  shell  to  generate  a 
set  of  velocity  response  curves.   The  velocity-time  history 
response  at  the  attachment  point  could  then  be  input  into 
the  SAPIV  code  for  the  finite  element  model  of  any  number  of 
small  pieces  of  equipment,  and  the  substructure  transient 
response  could  be  obtained  without  having  to  run  a  separate 
ELSHOK  calculation  for  each  shell/substructure  system.   This 
would  result  in  considerable  savings  in  computer  time. 


Although  ELSHOK  is  a  powerful  tool  which  can  be  used  to 
accurately  predict  the  results  of  underwater  explosion  tests 
on  equipment,  some  shortfalls  exist  in  the  program  which 
limit  its  use  in  design  applications.   At  the  present  time, 
there  is  no  direct  method  to  take  the  effects  of  the  surface 
cut-off  wave  or  bottom  reflection  into  account.   The  surface 
cut-off  wave  forms  in  shallow  underwater  explosions  when  the 
compressive  shock  wave  reflects  off  of  the  surface  of  the 
water  causing  a  tensile  wave  to  propagate.   Bottom 
reflections  form  in  deep  explosions  as  the  shock  wave  is 
reflected  off  of  the  floor  of  the  ocean.   ELSHOK  assumes 
that  the  explosion  is  occurring  in  an  infinite  acoustic 
medium  with  a  single  shock  wave.   Another  weakness  in  the 
program  is  the  fact  that  it  cannot  process  oblique 
explosions;  charges  must  be  placed  in  a  pure  end-on  or 
side-on  configuration.   Shock  tests  conducted  in  accordance 
with  MIL-S-901D  are  either  end-on  or  side-on  loading, 
however,  so  this  is  not  a  serious  drawback  in  predicting  the 
transient  motion  of  equipment  in  SSTV's. 

It  is  recommended  that  further  study  be  conducted  in 
this  area  using  more  refined  shell  models  resembling  full- 
scale  submarines.   Once  a  shell  model  is  produced  for  a 
given  submarine  class,  researchers  can  create  different 
finite  element  models  of  actual  equipment  to  observe  the 
interactions  between  the  submarine  hull  and  equipment  of 
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varying  mass  and  shape  when  subjected  to  end-on  or  side-on 
loads.   It  is  felt  that  ELSHOK  is  an  effective  method  to 
give  an  indication  of  the  response  of  equipment  in 
underwater  explosion  tests  without  needless  destruction  from 
unexpected  failures.   Its  continued  use  prior  to  UNDEX  tests 
will  lead  to  savings  in  both  time  and  effort. 
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APPENDIX  A 

PROGRAM  USED  TO  CONVERT  PUSLOB  PUNCH  CARD  FILES 
TO  FORMAT  ACCEPTANCE  BY  EASYPLOT 


Along  with  the  TEKTRONIX  plots  of  the  velocity-time 
histories,  PUSLOB  generates  velocity  punch  card  files  for 
each  requested  node.   The  punch  card  file  consists  of  a 
block  of  velocities  in  exponential  format  the  size  of  which 
depends  upon  the  number  of  time  points  requested.   The 
EASYPLOT  program  developed  for  the  IBM  3033  by  Mr.  John 
Mainwaring  is  a  quick  interactive  method  to  produce  graphs 
using  DISSPLA  without  having  to  create  a  plotting  program. 
It  will  accept  up  to  100  data  points  in  tabular  form 
provided  it  is  in  a  column  format  for  X-Y  input.   The 
PLOTCONVl  code  was  created  to  accept  data  from  two  separate 
punch  card  files  for  taper  and  conventional  charge 
comparisons  and  combine  them  into  one  data  file  for  each 
node  for  input  into  EASYPLOT. 

In  order  to  prevent  data  overflow,  the  time  step  must  be 
sized  in  USLOB  and  PUSLOB  to  generate  no  more  than  100 
velocity  data  points.   Both  the  PLOTCONVl  VAX  11/780  VMS 
computer  command  file  and  the  FORTRAN  code  are  included  in 
this  appendix.   The  command  file  must  be  modified  prior  to 
each  run  to  specify  the  punch  card  files  being  read  and  the 
titles  of  the  new  data  files. 
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$!  PL0TC0NV1.COM  -  -  COMMAND  FILE  FOR  EXECUTION  OF  PLOTCONVl 

$! 

$  ASSIGN  PUSL48.PUN  FOROll 

$    ASSIGN    PUSL«49.PUN    FOR012 

$  ASSIGN  MAXSHL1.DAT  F0R021 

$  ASSIGN  MAXSHL2.DAT  FOR022 

$  ASSIGN  MAXSHL3.DAT  FaR023 

$  ASSIGN  MAXSHL4.DAT  FOR02<+ 

$  ASSIGN  MAXSHL5.DAT  F0R025 

$  ASSIGN  MAXSHL6.DAT  FOR026 

$  ASSIGN  MAXSHL7.DAT  F0R027 

$  ASSIGN  MAXSHL8.DAT  FOR028 

$  ASSIGN  MAXSHL9.DAT  FOR029 

$  ASSIGN  MAXSHL10.DAT  FOR030 

$  ASSIGN  MAXSHL11.DAT  F0R031 

$  ASSIGN  MAXSHL12.DAT  F0R032 

$  ASSIGN  MAXSHL13.DAT  FOR033 

$  ASSIGN  MAXSHLl-i-DAT  FOR034 

$! 

$  RUN  PLOTCONVl 

$1 


PROGRAM  PLOTCONVl 
C 

C   THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  READ  TWO  PUNCH-CARD 
C   FILES  FROM  PUSLOB  AND  CONVERT  THEM  TO  A  FORMAT  WHICH  CAN 
C   BE  ACCEPTED  BY  EASYPLOT. 


C 


DIMENSION    VEL(100,1'4),VELO(100,3),VEL1(100,14),LABEL(14),TAG(14), 
%LTITLE(1<+) 

CHARACTERXIO  LABEL 

CHARACTER5«40  LTITLE 

COMMON/CSOLVE/DTRECS 

NCURV  =  1<+ 

DO  10  K=1,NCURV 

READ(11,999)  LTITLE(K) , NRECS, DTRECS 

READ(11,998)  LABEL(K),  TAG(K) 

READ(11,997)  (VELCJ,K),  J=I,NRECS) 

READ(12,999)  LTITLE(K) , NRECS, DTRECS 

READ(12,998)  LABEL(K) , TAG(K) 

READ(12,997)  (VEL1(J,K),  J=l, NRECS) 
10  CONTINUE 
999  FORMAT(A40,I5,1PE11.4) 
998  FORMAT(A10,F10.5) 
997  F0RMAT(1P7E11.4) 

DO  20  L=1,NCURV 

T  =  0.0 

LL=L+20 

DO  30  M=l, NRECS 

VEL0(M,1)=T 

DTRECS1=1000. OXDTRECS 

T=T+DTRECS1 

VEL0(M,2)=VEL(M,L)/12.0 

VEL0(M,3)=VEL1(M,L)/12.0 
30  CONTINUE 

DO  <i0  KK  =  1,  NRECS 

WRITE(LL,996)  (VELOCKK, J J ) , J J=l , 3) 
40  CONTINUE 

WRITE(LL,995)  LTITLEC L) , NRECS. DTRECS 

WRITE(LL,994)  LABEL ( L ), TAGC L ) 
20  CONTINUE 
996  F0RMAT(3F15.5) 
995  FORMAT(A40,I5,1PE11.4) 
994  F0RMAT(A10,F10.5) 

STOP 
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APPENDIX  B 
TYPICAL  ELSHOK  INPUT  CODES — SHELL  AND  SUBSTRUCTURE 

The  input  codes  which  follow  are  for  the  resonant 
diaphragm  case  for  both  taper  and  conventional  charges.   The 
only  differences  between  the  two  is  in  the  USLOB  code  where 
the  type  of  charge  is  specified.   All  the  input  codes  are 
explained  in  detail  in  the  ELSHOK  users  manual  [Ref.  4]. 

1.   B0S0R4  INPUT  DATA 

The  B0S0R4  input  code  is  created  by  working  inter- 
actively with  the  computer.   The  prompt  following  the  $ 
symbol  appears,  and  the  user  supplies  the  required 
information.   The  three  segments  (two  endplates  and 
cylinder)  are  modeled  separately,  and  then  they  are 
"connected"  in  the  global  data  section.   A  separate  input 
code  is  required  for  each  circumferential  harmonic  included 
in  the  calculation.   B0S0R4  output  provides  the  user  with 
the  in-vacuo  free-free  modes  and  natural  frequencies  of  the 
shell  . 
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FILE:  BOSSO 


DATA 


Al 


OSOR<*  ATTEMPT  ONE 


10 
3 
1 
0.0 
0.0 
16.31250 
0.0 
0 
3 
2.500000 
H 

0 
O.OOOQOOOE-t-00 
0 
0 
2 
0.1080000E+05 
0.3200000 
0.2535000E-03 
O.OOOOOOOE+00 
0 
1 
N 
N 
H 
H 
45 
3 
1 
16.31250 
O.OOQOOOOE-t-00 
16.31250 
253.1250 
0 
3 
0.1375000 
N 
6 
2 
3.500000 
53.12500 
103.7500 
15'i.3750 
205.0000 
254.6250 


0.3000000E+08 

7.125000 

7.535500 

2.375000 
Q.OOOOOOOE+00 
-1.781250 
O.OOOOOOOE+00 
0.7501E+a8 
0.7330000E-Q3 
0.3000000E+08 

3.250000 

11.69800 

2.750000 
O.OOOOOOOE+OO 


N  =  0 

INDIC  =  ANALYSIS  TYPE  INDICATOR 

NPRT  =  OUTPUT  OPTIONS  (1=MINIMUM,  2=MEDIUM,  3=MAXIMUM) 
ISTRES=  OUTPUT  CONTROL  ( 0=RESULTANTS .  1=SIGMA,  2=EPSIL0N) 
NSEG   =  NUMBER  OF  SHELL  SEGMENTS  (LESS  THAN  25) 

SEGMENT   NUMBER  11111111 

NMESH  =  NUMBER  OF  NODE  POINTS  (5  =  MIN.;  98  =  MAX.)C  1) 

NTYPEH=  CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 

NSHAPE=  INDICATOR  (1,2  OR  4)  FOR  GEOMETRY  OF  MERIDIAN 

Rl      =  RADIUS  AT  BEGINNING  OF  SEGMENT  (SEE  P.  66) 

21      =  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 

R2      =  RADIUS  AT  END  OF  SEGMENT 

Z2      =  AXIAL  COORDINATE  AT  END  OF  SEGMENT 

IMP    =  INDICATOR  FOR  IMPERFECTION  (0=NONE.  1=S0ME) 

NTYPEZ=  CONTROL  (1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 

ZVAL   =  DISTANCE  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

DO  YOU  WANT  TO  PRINT  OUT  R(S),  R'(S),  ETC.  FOR  THIS  SEGMENT? 

NRINGS=  NUMBER  (MAX=20)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 

IC=ELASTIC  FOUNDATION  MODULUS  (E.G.  LB/IN«X3)IN  THIS  SEG. 

LINTYP=  INDICATOR  (0,  1,  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 

NLTYPE=CONTROL  (0,1,2,3)  FOR  TYPE  OF  SURFACE  LOADING 

NWALL=INDEX  (1,  2,  4,  5,  6,  7,  3)  FOR  WALL  CONSTRUCTION 

E      =  YOUNG'S  MODULUS  FOR  SKIN 

U      =  POISSON'S  RATIO  FOR  SKIN 

SM  =MASS  DENSITY  OF  SKIN  (E.G.  ALUM.  =  .00025  LB-SEC3t3«2/INXK4) 

ALPHA  =  COEFFICIENT  OF  THERMAL  EXPANSION 

ANRS   =  CONTROL  (0  OR  1)  FOR  ADDITION  OF  SMEARED  STIFFENERS 

SUR    =  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

DO  YOU  WANT  TO  PRINT  OUT  THE  CCI.J)  AT  MERIDIONAL  STATIONS? 

DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 


SEGMENT 

NMESH  = 

NTYPEH= 

NSHAPE= 

Rl 

Zl 

R2 

Z2 

IMP    = 

NTYPEZ= 

ZVAL   = 

DO  YOU 

NRINGS= 

NTYPE  = 

Z(I)   = 

Z(I)   = 

Z(I)   = 

Z(I)   = 

Z(I)   = 

Z(I)   = 

NTYPER= 

NTYPER= 

NTYPER= 

NTYPER= 

NTYPER= 

NTYPER= 

E 

A      = 

lY 

IX 

IXY    = 

El 

E2 

GJ     = 

RM     = 

E 

A      = 

lY 

IX 

IXY    = 


NUMBER     22222222 
NUMBER  OF  NODE  POINTS  (5  =  MIN.;  98  =  MAX . ) (  2) 
CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 
INDICATOR  (1,2  OR  4)  FOR  GEOMETRY  OF  MERIDIAN 

=  RADIUS  AT  BEGINNING  OF  SEGMENT  (SEE  P.  66) 

=  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 

=  RADIUS  AT  END  OF  SEGMENT 

=  AXIAL  COORDINATE  AT  END  OF  SEGMENT 
INDICATOR  FOR  IMPERFECTION  (0=NONE,  1=S0ME) 
CONTROL  (1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 
DISTANCE  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

WANT  TO  PRINT  OUT  R(S),  R'(S),  ETC.  FOR  THIS  SEGMENT? 
NUMBER  (MAX=20)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 
CONTROL  FOR  IDENTIFICATION  OF  RING  LOCATION  (2=Z,  3=R) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  1) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  2) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  3) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  4) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  5) 
AXIAL  COORDINATE  OF  ITH  RING,  Z(  6) 


TYPE  (0  OR  1 

TYPE  (0  OR  1 

TYPE  (0  OR  1 

TYPE  (0  OR  1 

TYPE  (0  OR  1 


TYPE  (0 
YOUNG'S 


OR  1 


OR 
OR 
OR 
OR 
OR 
OR 


OR 
OR 
OR 
OR 
OR 
OR 


OR  5)  OF  DISCRETE  RING  NO . (  1) 

OR  5)  OF  DISCRETE  RING  N0.(  2) 

OR  5)  OF  DISCRETE  RING  NO . (  3) 

OR  5)  OF  DISCRETE  RING  NO . (  4) 

OR  5)  OF  DISCRETE  RING  N0.(  5) 


ON  P.70)(  1) 


70)(  1) 


OR  5)  OF  DISCRETE  RING  NO . (  6) 
MODULUS  OF  RING(  1) 
CROSS  SECTION  AREA  OF  RING(  1) 
MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG. 
MOMENT  OF  INERTIA  ABOUT  X-AXIS(  1) 
PRODUCT  OF  INERTIA(  1) 

RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P. 
AXIAL   COMPONENT  OF  RING  ECCENTRICITY(  1) 
TORSIONAL  RIGIDITY(  1) 

RING  MATERIAL  DENSITY  (E.G.  ALUMINUM= . 0002535) (  1) 
YOUNG'S  MODULUS  OF  RING(  2) 
CROSS  SECTION  AREA  OF  RING(  2) 

MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70)(  2) 
MOMENT  OF  INERTIA  ABOUT  X-AXIS(  2) 
PRODUCT  OF  INERTIA(  2) 
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FILE:  BOSSO     DATA     Al 

-2.062500      S  El  =  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70)(  2) 

(J.OOOOOQOE+QO  $  E2  =  AXIAL   COMPONENT  OF  RING  ECCENTRICITYC  2) 

0.92<*9000E+08  $  GJ  =  TORSIONAL  RIGIDITYC  2) 

0.7330000E-03  $  RM  =  RING  MATERIAL  DENSITY  (E.G.  ALUMINUM= . 0002535) (  2) 

0.3000000E+08  $  E  =  YOUNG'S  MODULUS  OF  RING(  3) 

■3.250000      «  A  =  CROSS  SECTION  AREA  OF  RING(  3) 

11.69800      $  lY  =  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70H  3) 

2.750000      «  IX  =  MOMENT  OF  INERTIA  ABOUT  X-AXIS(  3) 

O.OOOOOOOE+00  «  IXY  =  PRODUCT  OF  INERTIA(  3) 

-2.062500      $  El  =  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70)(  3) 

O.OOOOOOOE+00  $  E2  =  AXIAL   COMPONENT  OF  RING  ECCENTRICITYC  3) 

0.92<49000E+08  $  GJ  =  TORSIONAL  RIGIDITY(  3) 

0.7330000E-03  5  RM  =  RING  MATERIAL  DENSITY  (E.G.  ALUMINUM' . 0002535) (  3) 

0.3000000E+08  «  E  =  YOUNG'S  MODULUS  OF  RING(  <i) 

S. 250000      $  A  =  CROSS  SECTION  AREA  OF  RING(  ^i) 

11.69800      $  lY  =  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70)(  <i) 

2.750000      $  IX  =  MOMENT  OF  INERTIA  ABOUT  X-AXIS(  <i) 

O.OOOOOOOE+00  $  IXY  =  PRODUCT  OF  INERTIAC  <♦) 

-2.062500      $  El  =  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70)(  <♦) 

O.OOOOOOOE+OO  $  E2  =  AXIAL   COMPONENT  OF  RING  ECCENTRICITY(  4) 

0.92<i9000E+08  S  GJ  =  TORSIONAL  RIGIDITYC  ■*) 

0.7330000E-03  $  RM  =  RING  MATERIAL  DENSITY  (E.G.  ALUMINUM= . 0002535) (  <♦) 

0.3000000E+08  $  E  =  YOUNG'S  MODULUS  OF  RING(  5) 

8.250000      $  A  =  CROSS  SECTION  AREA  OF  RINGC  5) 

11.69800      «  lY  =  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  CSEE  FIG.  ON  P.70)C  5) 

2.750000      $  IX  =  MOMENT  OF  INERTIA  ABOUT  X-AXISC  5) 

O.OOOOOOOE+OO  $  IXY  =  PRODUCT  OF  INERTIAC  5) 

-2.062500      $  El  =  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70)C  5) 

O.OOOOOOOE+00  $  E2  =  AXIAL   COMPONENT  OF  RING  ECCENTRICITYC  5) 

0.92<i9000E+08  S  GJ  =  TORSIONAL  RIGIDITYC  5) 

0.7330000E-03  $  RM  =  RING  MATERIAL  DENSITY  (E.G.  ALUMINUM= . 0002535) C  5) 

0.3000000E+08  «  E  =  YOUNG'S  MODULUS  OF  RINGC  6) 

7.125000      $  A  =  CROSS  SECTION  AREA  OF  RINGC  6) 

7.535500      «  lY  =  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70)C  6) 

2.375000      *  IX  s  MOMENT  OF  INERTIA  ABOUT  X-AXISC  6) 

O.OOOOOOOE+OO  «  IXY  =  PRODUCT  OF  INERTIAC  6) 

-1.781250      «  El  =  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70)C  6) 

O.OOOOOOOE+00  «  E2  =  AXIAL   COMPONENT  OF  RING  ECCENTRICITYC  6) 

0.7501E+08     $  GJ  =  TORSIONAL  RIGIDITYC  6) 

0.7330000E-03  S  RM  =  RING  MATERIAL  DENSITY  (E.G.  ALUMINUM= . 0002535) C  6) 
O.OOOOOOOE+OO  $  K=ELASTIC  FOUNDATION  MODULUS  (E.G.  LB/INX*3)IN  THIS  SEG. 
0       $  LINTYP=  INDICATOR  (0,  1,  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 

0  S  NLTYPE=CONTROL  (0,1,2,3)  FOR  TYPE  OF  SURFACE  LOADING 

2  «  NWALL=INDEX  (1.  2,  ^,    5,  6,  7,  3)  FOR  WALL  CONSTRUCTION 
0.3000000E+08  $  E      =  YOUNG'S  MODULUS  FOR  SKIN 

0.3000000      S  U      =  POISSON'S  RATIO  FOR  SKIN 

0.7330000E-03  $  SM  =MASS  DENSITY  OF  SKIN  (E.G.  ALUM.  =  .00025  LB-SECX3i2/INXX<i) 
0.0       $  ALPHA  =  COEFFICIENT  OF  THERMAL  EXPANSION 

1  $  ANRS   =  CONTROL  (0  OR  1)  FOR  ADDITION  OF  SMEARED  STIFFENERS 
1       $  SUR    =  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

N  «  ARE  THERE  STRINGERS  (PLEASE  ANSWER  Y  OR  N)? 

Y  «  ARE  THERE  RINGS  (PLEASE  ANSWER  Y  OR  N)? 

0  S  K2   =CONTROL  (0  OR  1)  FOR  INTERNAL  OR  EXTERNAL  RINGS 
0.3000000E+08  $  E2     =  RING  MODULUS 

0. 3000000      $  U2     =  RING  POISSON  RATIO 
0.7330000E-03  $  RGMD=  RING  MASS  DENSITY 

Y  $  IS  THE  RING  CROSS  SECTION  CONSTANT  IN  THIS  SEGMENT? 

Y  $  IS  THE  RING  CROSS  SECTION  RECTANGULAR  (Y  OR  N)T 
5.625000      $  D2    =  ARC  LENGTH  BETWEEN  ADJACENT  RINGS  (CONSTANT) 

0.2500000  $  T2    =  THICKNESS  OF  RING  (CONSTANT) 

1.250000  $  H2    =  HEIGHT  OF  RING  (CONSTANT) 

N  $  DO  YOU  WANT  TO  PRINT  OUT  THE  CCI.J)  AT  MERIDIONAL  STATIONS? 

N  $  DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 

H  S 

H  $  SEGMENT  NUMBER     33333333 

10  $  NMESH  =  NUMBER  OF  NODE  POINTS  (5  =  MIN.;  98  =  MAX.)(  3) 

3  $  NTYPEH=  CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 

1  S  NSHAPE=  INDICATOR  (1,2  OR  4)  FOR  GEOMETRY  OF  MERIDIAN 
16.31250      $  Rl      =  RADIUS  AT  BEGINNING  OF  SEGMENT  (SEE  P.  66) 

0.0       «  21      =  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 
O.OOOQOOOE+00  $  R2      =  RADIUS  AT  END  OF  SEGMENT 

0.0       $  22      =  AXIAL  COORDINATE  AT  END  OF  SEGMENT 
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FILE:  BOSSO  DATA     Al 

0  $  IMP    =  INDICATOR  FOR  IMPERFECTION  (0=NONE.  1=S0ME) 

3  $  NTYPEZ=  CONTROL  CI  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 

2.500000  *  ZVAL   =  DISTANCE  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

N  9  DO  YOU  WANT  TO  PRINT  OUT  RCS),  R'CS),  ETC.  FOR  THIS  SEGMENT? 

0  $  NRINGS=  NUMBER  (:MAX  =  2a)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 

0.0  $  K=ELASTIC  FOUNDATION  MODULUS  (E.G.  LB/IN3«3<3)  IN  THIS  SEG. 

0  $  LINTYP=  INDICATOR  (0,  1.  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 

0  $  NLTYPE=CONTROL  (0,1,2,3)  FOR  TYPE  OF  SURFACE  LOADING 

2  $  NWALL=INDEX  (1,  2,  ^,  5,  6.  7,  3)  FOR  WALL  CONSTRUCTION 
0.1080000E+08  $  E      =  YOUNG'S  MODULUS  FOR  SKIN 

0.3200000  $  U      =  POISSON'S  RATIO  FOR  SKIN 

0.2535000E-03  $  SM  =MASS  DENSITY  OF  SKIN  (E.G.  ALUM. =.00025  LB-SECXX2/INKX<i) 

0.0  S  ALPHA  =  COEFFICIENT  OF  THERMAL  EXPANSION 

0  $  ANRS   =  CONTROL  (0  OR  1)  FOR  ADDITION  OF  SMEARED  STIFFENERS 

1  $  SUR    =  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -I) 

H  $  DO  YOU  WANT  TO  PRINT  OUT  THE  C(I,J)  AT  MERIDIONAL  STATIONS? 

N  $  DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 

H  $ 

H  $  GLOBAL  DATA 

-1  $  NLAST  =  PLOT  OPTIONS   (-1=N0NE,  0=GEOMETRY,  1=U,V,W) 

0  $  NOB    =  STARTING  NUMBER  OF  CIRC.  WAVES  (BUCKLING  ANALYSIS) 

0  $  NMINB  =  MINIMUM  NUMBER  OF  CIRC.  WAVES   (BUCKLING  ANALYSIS) 

0  $  NMAXB  =  MAXIMUM  NUMBER  OF  CIRC.  WAVES   (BUCKLING  ANALYSIS) 

1  $  INCRB  =  INCREMENT  IN  NUMBER  OF  CIRC.  WAVES  (BUCKLING) 
30  $  NVEC   =  NUMBER  OF  EIGENVALUES  FOR  EACH  WAVE  NUMBER 

O.OOOOOOOE+00  5  P       =  PRESSURE  OR  SURFACE  TRACTION  MULTIPLIER 
O.OQOQOOOE+00  $  TEMP    =  TEMPERATURE  RISE  MULTIPLIER 

O.OOOOOOOE+00  $  OMEGA   =  ANGULAR  VEL.  ABOUT  AXIS  OF  REVOLUTION  (RAD/SEC) 

I  $  NUMBER  OF  POLES  (PLACES  WHERE  R=0)  IN  SEGMENTC  1) 

H  $ 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.     I     1     I     1 

1  $  IPOLE  =  NODAL  POINT  NUMBER  OF  POLE,  IPOLE(  1) 

0  $  AT  HOW  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

N  $  IS  THIS  SEGMENT  JOINED  TO  ANY  LOWER-NUMBERED  SEGMENTS? 

0  $  NUMBER  OF  POLES  (PLACES  WHERE  R=0)  IN  SEGMENTC  2) 

H  $ 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.     2     2     2     2 

0  $  AT  HOW  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

Y  $  IS  THIS  SEGMENT  JOINED  TO  ANY  LOWER-NUMBERED  SEGMENTS? 

1  $  AT  HOW  MAY  STATIONS  IS  THIS  SEGMENT  JOINED  TO  PREVIOUS  SEGS.? 
1  $  INODE  =  NODE  IN  CURRENT  SEGMENT  (ISEG)  OF  JUNCTION,  INODEC  1) 
1  $  JSEG   =  SEGMENT  NO.  OF  PREVIOUS  SEGMENT  INVOLVED  IN  JUNCTION 

10  $  JNODE  =  NODE  IN  PREVIOUS  SEGMNT  CJSEG)  OF  JUNCTION 

1  $  IUSTAR=  AXIAL  DISPLACEMENT  C0=NOT  SLAVED,  1=SLAVED) 

1  $  IVSTAR=  CIRCUMFERENTIAL  DISPLACEMENT  C0=NOT  SLAVED.  1=SLAVED) 

I  $  IWSTAR=  RADIAL  DISPLACEMENT  C0=NOT  SLAVED,  1=SLAVED) 

1  «  ICHI   =  MERIDIONAL  ROTATION  C0=NOT  SLAVED,  I=SLAVED) 

O.OOOOOOOE+00  «  Dl     =  RADIAL  COMPONENT  OF  JUNCTURE  GAP 

O.OOOOOOOE+OO  $  D2     =  AXIAL  COMPONENT  OF  JUNCTURE  GAP 

Y  $  IS  THIS  CONSTRAINT  THE  SAME  FOR  BOTH  PREBUCKLING  AND  BUCKLING? 
1  $  NUMBER  OF  POLES  CPLACES  WHERE  R=0)  IN  SEGMENTC  3) 

H  $ 

H  $  CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.     3     3     3     3 

10  5  IPOLE  =  NODAL  POINT  NUMBER  OF  POLE,  IPOLEC  1) 

0  $  AT  HOW  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

Y  $  IS  THIS  SEGMENT  JOINED  TO  ANY  LOWER-NUMBERED  SEGMENTS? 

1  $  AT  HOW  MAY  STATIONS  IS  THIS  SEGMENT  JOINED  TO  PREVIOUS  SEGS.? 

1  $  INODE  =  NODE  IN  CURRENT  SEGMENT  (ISEG)  OF  JUNCTION,  INODE(  1) 

2  $  JSEG   =  SEGMENT  NO.  OF  PREVIOUS  SEGMENT  INVOLVED  IN  JUNCTION 
«5  $  JNODE  =  NODE  IN  PREVIOUS  SEGMNT  (JSEG)  OF  JUNCTION 

1  $  IUSTAR=  AXIAL  DISPLACEMENT  (0=NOT  SLAVED,  1=SLAVED) 

1  $  IVSTAR=  CIRCUMFERENTIAL  DISPLACEMENT  (0=NOT  SLAVED,  1=SLAVED) 

1  S  IWSTAR=  RADIAL  DISPLACEMENT  (0=NOT  SLAVED,  1=SLAVED) 

1  $  ICHI   =  MERIDIONAL  ROTATION  (Q=NOT  SLAVED,  1=SLAVED) 

O.OOOOOOOE+00  5  Dl     =  RADIAL  COMPONENT  OF  JUNCTURE  GAP 

O.OOOOOOOE+OO  $  D2     =  AXIAL  COMPONENT  OF  JUNCTURE  GAP 

Y  $  IS  THIS  CONSTRAINT  THE  SAME  FOR  BOTH  PREBUCKLING  AND  BUCKLING? 

Y  5  DO  YOU  WANT  TO  LIST  OUTPUT  FOR  SEGMENTC  1) 

Y  $  DO  YOU  WANT  TO  LIST  OUTPUT  FOR  SEGMENTC  2) 

Y  $  DO  YOU  WANT  TO  LIST  OUTPUT  FOR  SEGMENTC  3) 

N  $  DO  YOU  WANT  TO  LIST  PREBUCKLING  RESULTANTS  AND  RING  FORCES? 
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2.       ACESNID    INPUT    DATA 

The   ACESNID   computer    code   determines    the   virtual    mass 
array,    a   quantity   producing    the    late-time   contribution   of 
the    DAA.       Only    one    ACESNID   calculation    is    required    for    all 
circumferential    harmonics    of    the    shell.      The   bands   dividing 
the    areas    of    the    shell    are    identified    in    the    NLEFT,    NCYL, 
and    NRITE    entries,    while    the    fluid    properties    are    specified 
in    the    RHOFL   and   VSOUND    row. 


FILE:    ACEl  DATA  Al 

VIRTUAL    MASS    FOR    RING-STIFFENED   CYLINDER   WITH    FLAT    ENDS,    N=0,1 

0  1110/    NSTART,NFINIS,NFREQ,NVMASS,NCHECK 
15    51    15    3    0    /    NLEFT,NCYL, NRITE, NWBOSCNSYMF 

1  6    10    0    0    /    NORDER,NFENDS,NFCENT,NFCMPT,NOMIT 
9.5968'iE-5    5.833E4    0.005   /    RHOFL  ,  VSOUND,  ERR 
1111110111/    OUTPUT    FLAGS 

1.0   /    CPS(l) 

1  1  1  12  /  SEFS,JSGBEG,JPTBEG,JSGEND,JPTEND,LEFT 

2  1  2  47  /  SEFS,JSGBEG,JPTBEG,JSGEND,JPTEND, CENTRAL 

3  1  3  12  /  SEFS,JSGBEG,JPTBEG,JSGEND,JPTEND, RIGHT 

1  1  1  12  /  SOURCES, JSGBEG,JPTBEG,JSGEND,JPTEND, LEFT 

2  1  2  47  /  SOURCES, JSGBEG,JPTBEG,JSGEND,JPTEND, CENTRAL 
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3.       PIFLASH    INPUT    DATA 

The    PIFLASH   code    takes    data    from    the    shell    mode    files 
from  B0S0R4    and    the   virtual   mass    file    from   ACESNID   and 
combines    them    to    create    a    "shell-fluid    file"    containing    all 
the    information   concerning    the    shell    and    its    acoustic  medium 
for    further    calculations.      Since   only    the   breathing   modes 
are   activated    in   end-on    loading,    all    torsional    modes    are 
dropped    from    the    shell    file.      This    is    accomplished    in    the 
NJUSE    rows   where    the   modes    of    Table    I    are    identified. 


FILE:    PIEl  DATA  Al 

2  0    0    0    0/    NUMBER, NTORSN,NPTM,NSYMS,NSYMP 

3  3    3    3/    (NWETSG(K),NUSESG(K),K=1,2) 
386  .4    /    GRAVITY 

SHELLOO    /    SMF,    N=0 

SHELLOl    /    SMF,    N=l 

19    2   /    (NJUSE(J),J=1,NITEMS) 

12    3/    (KORSG(K),K=1,NKORSG) 

2    4    6    7    9    12    13    1*+    19    20    21    22    23    25    26    27    28    29    30   /    JUSE,    N  =  0 

1    2   /    JUSE,    N=l 
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4.   SAPIV  INPUT  DATA 

The  substructure  is  modeled  using  the  standard  SAPIV 
code  with  slight  modifications  by  the  developers  of  ELSHOK. 
The  diaphragm  is  modeled  using  twenty-four  plates  and 
twenty-five  nodes.   The  location  of  each  node  is  identified 
in  cartesian  coordinates.   Each  plate  is  1.06  in  thick  for 
the  resonant  case.   Nodes  eighteen  through  twenty-five  are 
rigidly  attached  to  the  shell  in  the  last  eight  lines  of  the 
code.   The  SAPIV  code  determines  the  fixed-base  modes  and 
corresponding  natural  frequencies  as  well  as  the  uncon- 
strained mass  and  stiffness  matrices  for  the  substructure. 
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FILE:  SAP23 


DATA 


Al 


IDEALIZED  INTERNAL  EQUIPMENT  AT  FRAME  9 


25  1  0  10 


0  1 
0  0 
0  0 


1 
2 
3 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
0  20  0 
0  21  0 
0  22  0 
0  23  0 
0  24  0  0 
0  25  0  0 


0  0  0  0 

1  1  0.0 
0  0  0.0 
0  0  0.0 

0  0.0 

0  0.0 


0 

0 

0 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0  0 
0  0 
0  0 


,N10SV,NRIGID 


/  NUMNP,NELTYP, 
0  0.0  0  0.0 
0  5.604  0  0.0 
963  3.963  0  0.0 
604  0.0  0  0.0 
963  -3.963  0  0.0  . 
0  -5.604  0  0.0 
-3.963  -3.963  0  0.0 
-5.604  0.0  0  0.0 
-3.963  3.963  0  0.0 
0.0  11.208  0  0.0 
7.925  7.925  0  0.0 
11.208  0.0  0  0.0 
7.925  -7.925  0  0.0 
0.0  -11.208  0  0.0 
-7.925  -7.925  0  0.0 
-11.208  0.0  0  0.0 
-7.925  7.925  0  0.0 
0.0  16.813  0  0.0 
11.888  11.888  0  0.0 
16.813  0.0  0  0.0 
11.888  -11.888  0  0.0 
0.0  -16.813  0  0.0 
-11.888  -11.888  0  0.0 
-16.813  0.0  0  0.0 
0  -11.888  11.888  0  0.0 
6  24  1  /  TYPE  6  (PLATE  AND  SHELL  ELEMENTS) 
1  .733E-03  0.0  0.0  0.0 

32.967E+06  9.89E+06  0.0  32.967E+06  0.0  11.538E+06 
0.0  0.0  0.0  0.0 


0  10 


0  0 
0  0 
0 
0 
0 
0 
0 
0 
0 


0  0  1 


0 

0 

0 

0 

8 

9 

2 

3 

4 

5 

6 

7  8  0  0 

10  11  3 

15  8  16  17  9  0  1 

16  9  17  10  2  0  1 

17  10  18  19  11  0 

23  16  24  25  17  0 

24  17  25  18  10  0 
0  0  0 
0.0  0 


/  ELEMENT  LOAD  MULTIPLIERS 


06 
06 


0  0 
0  0 


06  0  0 

06  0  0 


06  0  0 
06  0  0 
06  0  0  0 
06  0  0  0 
1.06  0  0  0 
0  1.06  0  0  0 

0  1.06  0  0  0 

1  1  1.06  0  0  0 
1  0  1.06  0  0  0 
1  0  1  .06  0  0  0 

0  0.0  0.0  0.0  0.0  0.0  /  CONC  MASSES 
0  0.0  0.0  /  ELEMENT  LOAD  MULTIPLIERS 


0  0  0  0.0  0.0  0  0.0  /  IFPR,IFSS, 


,NFO, SHIFT 


8  0  /  N 

18  1 

19  1 
1 
1 


20 
21 

22  1 

23  1 

24  1 

25  1 
END  0 


ACE,NSYMTY 


0  0 


DATA  (DUMMY  TITLE) 
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5.   PICRUST  INPUT  DATA 

The  PICRUST  code  takes  data  from  the  substructure  mode 
file  from  SAPIV  and  reorganizes  it  to  facilitate  the 
solution  of  the  transient  response  equations.   The  location 
of  the  substructure  attachment  points  to  the  shell  is 
identified  for  each  connecting  node. 


FILE:  PICRUSTl  DATA     Al 


1111011100 

0111110101/    OUTPUT    FLAGS    1-20 
2    0    0/    NJUSE,NHWSOB,NHWBAR 

18  2  11  0.0 

19  2  11  «+5.0 

20  2  11  90.0 

21  2  11  135.0 

22  2  11  180.0 

23  2  11  225.0 
2^  2  11  270.0 

25  2  11  315.0  /  NIPSUB,LBOSEG,LBOSPT,ANGDEG 
0.0  /  DEGROT 
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6.   USLOB  INPUT  DATA 

The  two  USLOB  input  codes  for  the  taper  and  conventional 
charges  appear  on  the  following  page.   USL46  DATA  is  the 
code  for  the  taper  charge.   It  contains  nine  discrete  data 
points  from  the  pressure-time  history  in  seconds  and  pounds 
per  square  inch.   There  are  1600  time  points  used  in  the 
numerical  integration  of  the  governing  equations  with 
information  saved  every  sixteenth  point,  so  only  100 
velocity-time  data  points  are  retained  for  EASYPLOT  display. 
The  charge  is  placed  along  the  negative  X-axis  at  840 
inches.   USL47  DATA  contains  the  information  for  the 
conventional  charge.   The  only  difference  between  the  two 
codes  is  a  charge  identification  entry  (NCHRG)  and  the 
charge  weight  and  constants  specified  for  the  conventional 
charge  instead  of  the  pressure-time  history  (WCHRG-THEXP) . 
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FILE:  USL^6 


DATA 


Al 


1599  16  2  9  2  2  1  /  NTIME, NSKIP, NCHRG, NQUAD, NFINE, KOUPLE, NSUBS 


OE-05  -840 

/  NSHAPE 

0  240.0  0.0002  320.0 

0034  115.0  0.0044  70 

0  /  DECAY 

10  0  0  0  0 

0  10  0  0  0 
/  NPTSHL 

1  2  11  2  46 
/  NPTSUB 


0  0.0  0.0  /  DELT,XLOAD,RLOAD,SURCUT 


00053  350 
0.0054  48 


00243  305 
00565  0.0 


0  0.00294  160.0 

/  (TC(J),PC(J),J=1, NSHAPE) 


0  0 

0  0 


/  OUTPUT  FLAGS  1-20 


3  12  /  LBOSEG,LBOSPT 


FILE:  USL47 


DATA 


Al 


1599  16  1  9  2  2  1  /  NTIME, NSKIP, NCHRG. NQUAD, NFINE, KOUPLE, NSUBS 

5.0E-05  -840.0  0.0  0.0  /  DELT,XLOAD, RLOAD, SURCUT 

352.0  3.8354314E+5  1.144  3.03131E-5  -.247  /  WCHRG, PZMLT, PZEXP,THMLT,THEXP 

0100000000 

0010000000/    OUTPUT    FLAGS    1-20 

4   /    NPTSHL 

1    1    2    11    2    46 

1   /    NPTSUB 


3    12   /    LBOSEG,LBOSPT 
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7.   PUSLOB  INPUT  DATA 

The  PUSLOB  code  produces  TEKTRONIX  plots  of  the 
velocity-time  histories  for  specified  points  on  the  shell  or 
substructure.   It  also  generates  velocity  punch  card  files 
which  are  converted  to  a  format  for  EASYPLOT  using  the 
PLOTCONVl  code. 


106 


FILE:  PUSL46    DATA     Al 


1599  16  1  1  1  /  NTIME,NSKIP,NSUBS,NTEK,NCARD 
5.0E-05  1000.0  1.0  /  DELT,XMULT,YMULT 

VELOCITY  PROFILE 

TIME  (MSEC) 

VELOCITY  (IN/SEC) 
2     NPTSHL-W 
V-FWD  1/01 

11-11    0.0 
V-FWD    3/12 

3    12    -1    1    0.0    /    LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 
13         NPTSHL-U 
V-FWD    1/01 

1  1    1    1    0.0   • 
V-FWD    2/02 

2  2    0    1    0.0 
V-FWD    2/07 

2   7    0    1    0.0 
V-FWD    2/11 

2    11    0    1    0.0 
V-FWD    2/15 
•      2    15    0    1    0.0 
V-FWD   2/20 

2    20    0    1    0.0 
V-FWD    2/2^ 

2    2<^    0    1    0.0 
V-FWD    2/28 

2    28    0    1    0.0 
V-FWD    2/33 

2    33    0    1    0.0 
V-FWD    2/37 

2    37    0    1    0.0 
V-FWD    2/^+1 

2    1^1    0    1    0.0 
V-FWD   2/46 

2  46    0    1    0.0 
V-FWD    3/12 

3  12    1    1    0.0    /    LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 
2  NPTSHL-V 

V-FWD    1/01 

11-11    0.0 
V-FWD    3/12 

3  12  -1  1  0.0  /  LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 
1      NPTSUB 
PLT  CNTR 
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